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Closed System Analysis 
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Abstract. Laser induced line narrowing effect, dis- 
covered more than thirty years ago, can also be applied 
to recent studies in high resolution spectroscopy based 
on electromagnetically induced transparency. In this pa- 
per we first present a general form of the transmission 
width of electromagnetically induced transparency in a 
homogeneously broadened medium. We then analyze a 
Doppler broadened medium by using a Lorentzian func- 
tion as the atomic velocity distribution. The dependence 
of the transmission linewidth on the driving field inten- 
sity is discussed and compared to the laser induced line 
narrowing effect. This dependence can be characterized 
by a parameter which can be regarded as "the degree of 
optical pumping" . 

PACS: 32.70.Jz, 42.50.Gy, 42.55.-f, 42.65.-k 



Over the last decade, considerable attention has been 
paid to the studies of the atomic coherence effects and 
their applications Eip. The technique of Electromag- 
netically Induced Transparency (EIT) which makes an 
opaque medium become transparent by applying an ex- 
ternal coherent radiation field yields various appli- 
cations from enhancement of nonlinear optical processes 
to slow light H^Q. In addition to the elimina- 
tion of absorption, the absorption profile reveals a narrow 
transmission line, which has been applied to high reso- 
lution spectroscopy and high sensitivity magnetometer 

Since many of these experiments are performed in an 
atomic cell configuration, the Doppler broadening effect 
on EIT is an important concern. Recent theoretical inves- 
tigations of Doppler broadening effects on EIT, however, 
has been focused mainly on the existence of EIT for cer- 
tain configurations Jl9|-pl|. The issue of EIT linewidth 
for a Doppler broadened medium has been lately ad- 
dressed by Taichenachev and coworkers |22] |. As the 
width of transmission line is directly related to the dis- 
persion near the EIT resonance, it is also a key issue in 
dispersive measurements. 



In a three-level A-type system if the system is homoge- 
neously broadened, as is well known, EIT can be achieved 
when the intensity of the driving field {fl 2 ) is larger than 
the product of the decay rate of the coherence between 
the lower levels ('Jbc) and the homogeneous linewidth (7). 
Then, if the system is inhomogeneously broadened (say, 
with the width Wd), one might guess that EIT can be 
achieved when fl 2 is larger than ^cWd instead of jbcj- 
This is not so. We show that one can still have EIT when 
2> 76c7 even in the case of inhomogeneous broaden- 
ing. 

For the spectral width of EIT, if the system is homo- 
geneously broadened, the two absorption lines are sepa- 
rated approximately by the Rabi frequency of the driving 
field f2 when Q, is larger than the homogeneous linewidth 
7. When Q, <C 7, it becomes f2 2 /7- Then, if the sys- 
tem is inhomogeneously broadened, it might be inferred 
that the EIT width goes as £1 when Vl is larger than the 
inhomogeneous linewidth Wd, and becomes il 2 /Wd as 
ft< W D - 

In the literature, however, we find that the narrow fea- 
ture superimposed on the Doppler broadened profile has 
been studied more than thirty years ago. Laser induced 
line narrowing effect was discovered by Feld and Javan 
p3| and the spectral width of the narrow line was shown 
to be linearly proportional to the driving field Rabi fre- 
quency. Various aspects of this effect has been investi- 
gated by Hansch and Tosche k |24[ , and it was also called 
nonlinear interference effects J25[. In a recent article [ p6[ , 
it has been proposed that this laser induced line narrow- 
ing can be applied to the recent experiments based on 
EIT and the spectral line of the EIT resonance can be 
narrower in a Doppler broadened system than in a ho- 
mogeneously broadened system. Here we analyze these 
ideas in detail and demonstrate the power broadening of 
the linewidth of EIT resonance in a Doppler broadened 
system. 

Under the condition of fl <C Wd , there are again two 
different regimes of EIT width: In one limit it is propor- 
tional to the Rabi frequency of the driving field, which 
has the same expression as the spectral width shown in 
the study of laser induced line narrowing p3|. As the 
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driving field gets strong, it becomes power broadened 
and indeed has a form proportional to the intensity of 
driving field (as £1 2 /Wd)- 

This paper is organized as follows: In Sec. I we set 
up our model scheme of the three-level system and the 
transmission width of EIT in a homogeneously broadened 
medium is discussed. In Sec. II the Doppler averaged sus- 
ceptibility is obtained by using a Lorentzian function for 
the velocity distribution and the absorption profile, the 
EIT condition, and the linewidth of EIT are discussed. 
Comparison between the closed system and the open sys- 
tem is briefly given in Sec. III. Section IV contains the 
summary of the present paper. 



I. HOMOGENEOUSLY BROADENED SYSTEM 

We consider a model scheme depicted in Fig. 1. The 
transition a «-> c is coupled to a coherent driving field and 
the transition a <-> b is coupled to a weak probe field. The 
atom-field interaction Hamiltonian can be written as 



V 



-hae 



! |a)(fe|-^e- W0 »(c|+H.c 



(1) 



where a is the Rabi frequency of the probe field, fl are 
the Rabi frequency of the driving field. In this model 
we take the decay rate from the level a to b (c) as 7 
(7'). The relaxation between the lower levels is denoted 
as 7& c such that the decay rate of the off-diagonal density 
matrix element (pbc) is defined as 7 bc . 




FIG. 1. Three-level model scheme. The upper level a de- 
cays to b and c with decay rate 7. The relaxation rate between 
levels b and c is denoted as "/be which is assumed to be small 
compared with 7. 

The equations of motion for the density matrix ele- 
ments in a rotating frame are then given by 

pab = -^abPab - ia(p aa - Pbb) + i^Pcb (2a) 
Pcb = -FcbPcb - iapca + iQpab (2b) 

Pac = -^acPac - iapbc - iQ(Paa - Pec) (2c) 
Pec = -JbcPcc + IPaa + JbcPbb - iQ(p C a ~ Pac) (2d) 
Paa = -(7 + l')paa ~ ia{Pab ~ Pba) ~ «^(Pac ~ Pea)- (2e) 



-fab = lac = ~ (7 + 7' + 76c), 



7cb = Ibc- 



(3) 



and Ay's are given as A a6 = Lu ab - v, A QC = u ac - u , 
and A cb = A afc - A ac . 

For a weak probe field, first order solution for the off- 
diagonal density matrix element p a b (which governs the 
absorption of the probe field) can be found in steady state 
as 



P {1) 

Pab 



FabFcb + f^ 2 



r cb (p. 



(0) 
aa 



J- c 



(4) 



where p^ is the population in level I in the absence of 
the probe field. The susceptibility is then written as 



X = V 




(5) 



where 77 is given by 77 = (3/87r)A^7A 3 for the atomic num- 
ber density N and the wavelength A. The effect of the 
probe field intnsity on the susceptibility is ignored by 
using the linear approximation p2]. 



A. Optical pumping and population distribution 

Let us find the population of each level in the absence 
of the probe field (i.e. the zeroth order population). Ob- 
viously, if the driving field is not turned on, we have 
Paa = and pbb — Pec — 1/2 from Eq. (^). Now as the 
driving field being turned on, in steady state, we have 
from Eq. (||c,e) 



p (0) 
rac 



p (0) 

raa 



' -p (.Paa Pcc^)-! 

J- n 



■ ac 

in 



7 + 7' 



/ \rac rca J 



(oh 



(6) 



Let us now assume, for the sake of simplicity, that the 
decay rate from the level a to c is same as the decay 
rate from the level a to b, i,e. 7' = 7 and the driv- 
ing field detuning is denoted as Ao. Then, we have 
r QC = 7a C + «A QC = (27 + 7 &c )/2 + iA Q , and 



(0) _ (o) = -»2n(7 + 7be /2) / (0) _ (0) 

/ ac Pea ^ + 7fcc/2)2 + A 2 [Paa Pec 



By Eqs. (||) and ([?]) we obtain 

27 



Q2 
X 



o 2 

JO) = » (0) 

raa rcc ' 



(7) 



(8) 



where 



Here we assume that the Rabi frequencies are real, 's 
are defined as 7^ + lAy, where 



X 



[(7 + 76c/2) 2 + Ag] 
2( 7 + 7bc /2) 



(9) 
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Note that Eq. (||d) can be written as 

Pec = - \ 7bc + ~—^J Pec + ^7 + -j^J Paa + IbcPbb (10) 

Hence using p^b = 1 — p a a — Pec, we obtain the zeroth 
order population 



(0) _ 2-lbcQ 

Paa ~~ 



(0) 
Pbb 



2D 1 

4 7 X 7fcc + 2 7fee ft 2 + 2ft 2 7 
2D 



(0) _ ±lX~i bc + 27 fee » 2 



2D 



(11a) 
(lib) 
(11c) 



where D = A-jbclX + 37& c ft 2 + ft 2 7- For 7 ^> 7 & c , these 
can be simplified as 

(Q) n (Q) ~ 4 7^7bc (o) _ 4 7 X7 be + 2ft 2 7 

aa Pec ~ ' ^»£> ~ 2D ' 



where 



x M 7 2 + A 2 ^ £ « 4 7bc7 X + ft 2 7 . (13) 
27 

Note that when the driving field is on resonance, the 
usual EIT condition ft 2 3> 7fc c 7 is equivalent to pbb ~ 1 
in Eq. (|l2|); i.e. a complete optical pumping to the level 
b is required to achieve EIT. 



A 4 -A 2 (2ft 2 + 7 2 ) + ft 4 = 0, 
and the solution is 



A = — 

2 



2s + 1 ± V4s + 1 



(17) 



(18) 



where s = ft 2 / 7 2 . Hence for s > 1 we have 



A 



7 2 (s + v 7 ^) 



a w ± n± 



2' 



(19) 



which shows that the absorption peaks are at ±ft with 
full width 7 and the half width of transmission is obtained 
as 



EIT 



ft 



7 



On the other hand, for s <C 1 we have 
,,2 r 



7 



2.s + 1 ± (1 + 2s- 2s 



Therefore, 



A«± 7 



ft 2 

7 



±- 



n 2 

7 



(20) 



(21) 



(22) 



B. Transmission width of EIT 

Now let us consider the transmission width under the 
condition of a resonant driving field. When we have a 
resonant driving field, i.e. A = 0, from Eq. ( |l3) . we 
find X w 7 /2 and D » ft 2 7 . Therefore, p£° 5 « pi? ~ 
and p^/ « 1, i.e. all the populations are in the level 
b. As is discussed in the previous section, the condition 
ft 2 76c7 leads to a complete optical pumping in the 
homogeneously broadened case. 

Eqs. (IH then yield 



X 



7?H)r c &(-i) 



r a fcr c (, + ft 2 

Since T a ft 7 + iA and r c ;, = j^c + *A, we have 



m . . . 

X = ^(7bc + J A) 



(ft 2 



iA 7 



(14) 



(15) 



where Z = (ft 2 — A 2 ) + A 2 7 2 . Hence, the imaginary part 
is obtained as 



n V 

X =Z 



76c (ft 2 - A 2 ) 



A 2 7 



(16) 



Since the maximum of x" is l/ 7 at A « ft, we may define 
T EIT , the half width of EIT as x"( A = T EIT ) = l/2 7 , 
which gives 



Hence, when ft <C 7, we have the absorption profile show- 
ing a whole envelope with half width 7 + ft 2 /7, and at 
the center there exists a transmission line with its half 
width as 



r 



EIT 



ft 2 / 7 . 



(23) 



We note that under the EIT condition ft 2 ^> 776c, ^eit 
cannot be smaller than 7 f, c . 



II. INHOMOGENEOULY BROADENED SYSTEM 

Now if the system is Doppler broadened, for the atoms 
with velocity v, the radiation fields are Doppler shifted 
as v — > v(\ — v/c) = v — kv for the probe field with k 
as the component of the wavevector on the propagation 
axis, and vq — > uq{X — v/c) — i/q — k'v for the driving 
field. Hence, for a Doppler broadened system, we re- 
place Ay as A q6 — * A Q& -I- kv A ac — > A QC + k'v, and 
A c b — > A c (, + (k ~ k')v. In the present analysis we as- 
sume that the energy difference between the level b and c 
is small enough so that we have k' sa k and the probe field 
and the driving field are copropagating such that (fc— k )v 
term can be neglected. Hence the atomic polarization 
should be averaged over the entire velocity distribution 
such that 
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X = / d(kv) f(kv) 77 



Pabikv) 



(24) 



where f(kv) is the velocity distribution function, and 
again r\ is given by Eq. (|J). We now consider the case 
where the inhomogeneous line is bigger than any other 
quantities involved such that Wd S> fi, 7 3> jbc, and the 
condition Jl 2 ^> 7^7 is still satisfied. 

The population distribution in Eq. (|l^) is now dif- 
ferent for atoms with different velocities. As we men- 
tioned in Sec. II, we then need to replace A with 
Ao + k'v 1=3 Ao + kv for the expression of X in Eq. Jl3| ) 
such that for a resonant driving field (Aq = 0) we have 



X 



7 2 + (kv) 2 
2^ 



D^2 lbc [j 2 + (kv) 2 ] + n 2 1 . (25) 



Hence, for the atom with its velocity v, p a b(kv) can be 
written 



pab(kv) 



ia 1 
Y 2D 



r cfc (4 7 X 76c + 2fi 2 7 ) 

^ 2 4 7 ^7bc 
7 + 7fc c /2 + ikv 



(26) 



where Y = (7 + 7i, c /2 + iA + i/cv)(7 bc + zA) + J1 2 . Here we 
have assumed k' w k and (k — k')v terms can be neglected 
for the copropagating fields. 



A. Doppler average using a Lorentzian distribution 

We now need to evaluate the expression of susceptibil- 
ity given in Eq. (|24|). Normally, the velocity distribution 
is described by a Gaussian function given by 



f{kv) 



1 



y/nki 



■ exp 



(kv 



,21 



(ku) 



(27) 
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FIG. 2. Velocity distribution of FWHM 2W D = IOO7 as 
a function a kv in unit of 7, with (a) a Gaussian profile of 



Eq. 
Eq. 



27), (b) a Lorentzian profile of Eq.(|28|), (c) the plot ofr 
2£ ) multiplied by a factor \/ n In 2. 



In Fig. 3 the absorption profiles are described numerically 
by using the two different distributions. We note that the 
two distributions give an almost identical result when the 
factor \Jn In 2 is taken into account, see Fig. 3(c). 
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FIG. 3. Absorption profiles (x^/v) as a function of probe 
field detuning (A in unit of 7) for 2Wd = IOO7, 7b c = 10 _3 7, 
and SI = 27 using (a), (b), (c) of Fig. 2, respectively. 



where u = \j2ksT jM is the most probable speed of 
the atom given by temperature T and the atomic mass 
M. Then, the full width at half maximum is given as 
2Wd — 2Vhi 2ku. However, in our analysis, for the sake 
of simple analytic expressions, we adopt a Lorentzian dis- 
tribution of FWHM of 2Wd, instead of a Gaussian dis- 
tribution, such that 



f(kv) 



Wd/tt 



W 2 



(kv) 2 



(28) 



The two distributions are shown in Fig. 2, and there 
we see that a Gaussian distribution with the same width 
(FWHM 2Wd) has its maximum larger than that of 
Lorentzian distribution by a factor of \/ tt In 2. Hence, 
if we multiply the factor v^rTn2 in Eq. (^), the central 
distribution becomes very similar to that of Gaussian as 
illustrated in Fig. 2(c). 



Now using the distribution of Eq. (|28|), Eq. ( |24| ) may 
be considered as a contour integration in the complex 
plane. We find three poles in the the upper half plane at 



kv = 



A($l 2 - 7fc 2 c - A 2 ) + z(A 2 7 + 7fc 2 c 7 + 7^ 2 ) 




iW D , 



and two poles in the lower half plane as 



(29) 



iWr 



'Jl 2 7 
276c 



(30) 



We can see that one pole is from the expression Y, 
two poles (±«y / Jl 2 7/27f )c ) are from the expression D in 
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Eq. (|26|), and two poles (±iWo) are from velocity distri- 
bution function |27|] . Let us take the contour in the lower 
half plane and denote 

X = Xi+X2, (31) 

where x^'s are the contributions from the two poles at 
—iWo and —iy / fl 2 j/2j bc , respectively. For the pole at 
kv = —iWjj, we obtain 



Xi 



-irj 



{Bi 



2Z X A 
where A is given by 



iAW, 



D 



Ci - iADi 



and 



A = -2 lbc W 2 D + fi 2 7 , 



Zi = {jbcW D +fl 2 - A 2 ) 2 + A 2 Wl, 

B 1=lbc w D + n 2 , 

a = 2 lbc w D { lbc w D + n 2 ) - 2 lbc n 2 1: 

Dx = -2 lbc W 2 D + 2fi 2 7 . 



For the pole at kv = —i^/fl 2 j/2'y bc , we have 
ir/il 2 jWD 



X2 



2Z 2 Ay 



(B 2 



iAy 



C 2 -iA 



where y = \/0 2 7 /2"fbc, and 

Z 2 = {ibcV + « 2 - A 2 ) 2 + A 2 y 2 , 
B% = j bc y + fi 2 , 
C 2 = -7&c + n 2 /y- 



(32) 



(33) 



(34) 



(35) 



(36) 



Note that we have assumed fi 2 ^> 7b c 7, Wd ^> 12, 7 ^> 

Ibc- 



X"(A = 0) 



Vlbc 



lbc w D + n 2 



l + Vx 



(39) 



where 



n 2 1 

2lbcW 2 • 



(40) 



T\\fx ri , , 



We note that when x -C 1, 

and when i>1, 

"i ?77bc ?77bc ?77 / ^ q\ 

x U=o => + n i < « (42) 

In both cases the EIT can be achieved, i.e., x"\a=o 
rj/W/j. Therefore, the condition for EIT is still Q 2 » 
77b c , the same as in the homogeneously broadened sys- 
tem. 

One interesting quantity here is the slope of the real 
part of the susceptibility, which is important in preci- 
sion magnetometry, and also governs the group velocity 
of the probe light. From Eqs. ([32 35) the real part of the 
susceptibility is found as 



Xi 

X2 



-77A 



2ZxA 

r]n 2 -fW D 



WdC\ 4- D\{B\ — A 2 ) 



C 2 y + (B 2 -A 2 ) 



2Z 2 Ay 

and its derivative at resonance is given by 



(43) 



B. Absorption and dispersion at EIT resonance 

The absorption profile is now obtained by the imagi- 
nary parts of Eqs.(P4 B5h as 



Xi 



-v 



2Z X A 



x 2 - 



2Z 2 Ay 
Taking A = 0, we found 

V 



(Bi - A 2 )Ci - A 2 W D D! 
(B 2 - A 2 )C 2 - A 2 y 



X ' 1 '(A = 0) = -^ 



X2(A = 0) 



lbcW D 

•nibcWo 



Ibcti 2 ! 



A 



1 



VlbcW D + Q 2 
ZjbcV 



j bc y + 12 2 



(37) 



(38) 



which gives the minimum value of absorption at the EIT 
line center as 



OA 



A=0 



777 

" A ' 



dxk 
dA 



^2 lbcl w 2 D /n 2 



A=0 



.4 



Hence, we obtained the slope of x' & t A = as 



dA 



A=0 



J_ 

ft 2 1 



(44) 



(45) 



Therefore, when x 3> 1, it approaches to rj/il 2 and when 
x -C 1, it goes as {■q/^l 2 )^/x. We note that, under the 
EIT condition 12 2 ^> 7f, c7 , iv/^ 2 )^ is still much larger 
than (r;/0 2 )(7/VF D ). 



dA 



1 



A=0 



(46) 
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C. Transmission width of EIT resonance 

In order to estimate the linewidth of EIT we take the 
same procedure as in Sec. II: First, we find that the max- 
imum of x" as Xmax ~ f]/WD at A « ±0. Then, we 
evaluate A which defines Veit as 

X "(A = T EIT ) = rj/2W D . (47) 

By Eq. (^7|) it readily gives the following equation: 

A4 _ 2 lbc &2 lbc wl + n 2 j A2 _ 2 7bc » 2 ^ 4 w2 = Q 



7 



2 lbc W L 



7 



(48) 



which yields the half width of the EIT for the Doppler 
broadened system given by 



r 2 

1 EIT 



7 

27b, 

7 



fi 2 (l +x) 



n 2 (i + x), 



Ax 



{l + xf 



1/2' 



(49) 



where x — Q 2 ^/2^/ bc W D given by Eq. ([h)|). Now if we 
define a saturation intensity as 



2 lbc W 2 D 

7 



the linewidth expression can be written as 



27fee 

7 



n 




(50) 



(51) 



Here we can see that in the limit f2 <gC T E it is propor- 
tional to the Rabi frequency of the driving field. Such a 
linewidth was predicted by Feld and Javan in the study 
of laser induced line narrowing p3[ . On the other hand, 
in the limit il 3> Sl s Teit is proportional to the inten- 
sity of the driving field (SI 2 /Wd)- This power broadening 
feature is shown in Fig. 4. 
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A 

FIG. 4. Absorption profile for 2Wd = IOO7, 7 6c = 10~ 3 7. 
(a) Q = 7, (b) n = 37, (c) Q = 67. Note that y^Tf ~ 0.03 7 
and Sl s ~ 4.5 7 . 



The expression of Eq. (|51j) shows a reminiscence of 
power broadening factor in the description of hole burn- 
ing [^8| . In place of the homogeneous linewidth in the ex- 
pression of hole burning, here we have an effective width 
which is determined by the spectral packet involved in 
population trapping pt|. 



D. The role of optical pumping 



We have seen that the parameter x — Sl/f2 s plays an 
important role in the case of inhomogeneously broadened 
medium. Let us here examine the physical meaning of the 
parameter. 

Suppose the system is homogeneously broadened. 
When the driving is on resonance, the optical pumping 
rate from the level c is then order of f2 2 /7, as given in 
Eq. (|l0|). A complete optical pumping within the homo- 
geneous linewidth, is then possible if this rate is larger 
than the pumping from level b to c: f2 2 / 7 3> 7 & c - This, in 
turn, gives the EIT condition. When we have the driving 
field detuned by A , the optical pumping rate decreases 
by a factor of 7 2 /( 7 2 + A 2 ,). Again for a complete optical 
pumping we need Jl 2 7 /[7 2 + A 2 ] ^> 7 hc . 

If we now assume that we have the resonant driving 
field A = again, and, instead, the atoms are moving. 
Then, for atoms with velocity v, the optical pumping 
rate becomes Sl 2 7 /[7 2 + (kv) 2 ]. Then, on the average, 
to have a complete optical pumping in a Doppler broad- 
ened system we need to require Jl 2 7 /( 7 2 + W73) ^> 7 b c , 
which corresponds to 1 > 1 (assuming Wd ^S> 7), i.e. 
ft ^> Q s = 2^ bc W D /-f. Hence, the parameter x = Sl 2 /f7 2 
represents the degree of saturation in b <-» c transition, 
or the degree of optical pumping from the level c to b 
within the inhomogencous linewidth. 



III. COMPARISON WITH AN OPEN SYSTEM 
DESCRIPTION 



In this section we examine the case of an open system 
and show that the result is essentially the same as our 
model of a closed system. The open system is modeled 
for the atoms that are coming in and out of the interac- 
tion (with the radiation fields) region. Although in such 
a case all the levels have the same decay rate (say, 7& c ), 
the upper level can decay much faster than the time of 
flight through the interaction region (for example, radia- 
tive decay or collisional decay). Hence we assume that 
the lower levels b and c decay with rate ^ bc and the upper 
level a decays with rate j a which is much bigger than j bc 
(see Fig. 5). 



G 




b The 



FIG. 5. Model scheme of the open system. Upper level a 
decays with rate 7„ Lower levels b and c decay with the same 
rate 7(, c . Atoms are pumped at a rate r equally to the lower 
levels. 

Furthermore, for simplicity, we assume that the atoms 
are coming into the interaction region with a same rate 
for the lower levels. Under these assumption, the equa- 
tion of motion for the density matrix elements can be 
written as 



have the open system description almost identical to the 
description for our model scheme of the closed system 
A detailed analysis of the open system will be presented 
elsewhere ifEJl. 



IV. SUMMARY 

In this paper, we have studied the transmission width 
of EIT in a three-level A system. The Doppler averaged 
susceptibility is found by using a Lorentzian velocity dis- 
tribution rather than the Gaussian distribution. Then 
we have shown the requirement for achieving EIT, and 
the analytic expression of the EIT linewidth. The satu- 
ration intensity f2 2 defines the degree of optical pumping 
as fi 2 /f2 2 , and represents the condition under which the 
broadening is either linear or quadratic in the Rabi fre- 
quency of the driving field. 



Pab = -^abPab ~ iOi{p a a - Pbb) + i^Pcb, (52a) 

Pcb = -FcbPcb - ictpca + iQpab, (52b) 

pac = -Tacpac - ictpbc - i^{paa - Pec), (52c) 
Paa = ~laPaa ~ ia(p a b - Pba) - i^iPac - Pea), (52d) 
Pbb = r - "fbcPbb - ia(pab - Pba) - i$l{Pae - Pea), (52e) 

Pcc = r--f bc p cc -iil(p ca - p ac ). (52f) 

Here the notations are the same as Eq. (||). Note that 
now we have r ac = (~f a + jbc)/2 + *Aq, which gives 



X' 



[(7a/2 + 7ftc/2) 5 



A 2 ,] 



2( 7a /2 + 7bc /2) 



(53) 



Again if we assume that j a 3> "fbc, the populations are 
found as 



p(o) 

raa 



where 



IbelaX' 

2D 1 ' 



(o) 
Pbb 



IbelaX' + » 2 7a 

2D 1 



(54) 



X' 



(7a/2) 2 



A 2 



7a 



D' W JbelaX' + n 2 la . (55) 



Comparing Eq. (|i|) with (|T^), we can see that the pop- 
ulation distribution is almost identical to the one for the 
model of closed system. 

Furthermore, the expression for p£? is identical to the 
one for the closed system given in Eq. (Eh. Let us then 
recall Eq. (|33|) saying that A = -2-y bc wl + tt 2 -f, which 
is obtained by putting — iWu to Ao in the expression 
of D in Eq. (JlT|). The sign of A determines whether 
the crucial parameter x is > 1 or < 1. Similarly, here 
for the open system, when we put Ao = —iWo, we can 
define A' as A' w —jbcWp + f2 2 7a such that we have 
x 1 = Qja/jbcWp as the parameter which plays the same 
role as x in Eq. (Il3|). Hence, by replacing 7 a ^> 2 7 , we 
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